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Summary. - The Lagrangian formulation of classical field theories and in particular 
general relativity leads to a coordinate-free, fully covariant analysis of these constrained 
systems. This paper applies multisymplectic techniques to obtain the analysis of Palatini 
and self-dual gravity theories as constrained systems, which have been studied so far in 
the Hamiltonian formalism. The constraint equations are derived while paying attention 
to boundary terms, and the Hamiltonian constraint turns out to be linear in the mul- 
timomenta. The equivalence with Ashtekar's formalism is also established. The whole 
constraint analysis, however, remains covariant in that the multimomentum map is eval- 
uated on any spacelike hypersurface. This study is motivated by the non-perturbative 
quantization program of general relativity. 
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1. - Introduction. 

When Dirac developed his approach to constrained Hamiltonian systems with the 
corresponding quantization program, he emphasized that the Hamiltonian formalism is 
always necessary to quantize a field theory with constraints. Remarkably, he provided 
a well-defined (though not unique) prescription to define a Hamiltonian function on the 
whole phase space, a classification of constraints in terms of their Poisson brackets which 
is immediately relevant for quantization, and an approach to quantum electrodynamics 
which does not rely on Feynman diagrams and renormalization theory [1-5] . However, any 
attempt to combine Dirac's quantization of first-class constrained Hamiltonian systems 
with the Arnowitt-Deser-Misner geometro dynamical framework for canonical gravity faces 
very severe technical problems. In other words, the occurrence of the scalar curvature 
of the spacelike three-surfaces and the product of the three-momenta in the Hamiltonian 
constraint make it impossible to find exact solutions of the corresponding Wheeler-De Witt 
equation, as well as interpret the (as yet unknown) physical states of the quantum theory 
within the geometrodynamical framework [5,6]. 

More recently, the work by Ashtekar, Rovelli, Smolin and their collaborators on con- 
nection dynamics and loop variables has made it possible to cast the constraint equations of 
general relativity in polynomial form, and then find a large class of solutions to the quan- 
tum version of constraints [7-14]. However, the quantum theory via the Rovelli-Smolin 
transform still suffers from severe mathematical problems in 3+1 space-time dimensions 
[15], and there appear to be reasons for studying non-perturbative quantum gravity also 
from a Lagrangian, rather than Hamiltonian, point of view (see below). The aim of this 
paper is therefore to provide a multisymplectic, Lagrangian framework for general rel- 
ativity [16-18], to complement the present attempts to quantize general relativity in a 
non-perturbative way. The motivations of our analysis are as follows. 

(i) In the case of field theories, there is not a unique prescription for taking duals, on 
passing to the Hamiltonian formalism. For example, algebraic and topological duals are 
different. In turn, this may lead to inequivalent quantum theories. 
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(ii) The 3+1 split of the Lorentzian space-time geometry, with the corresponding E x R 
topology, appears to violate the manifestly covariant nature of general relativity, as well 
as rely on a very restrictive assumption on the topology [19]. 

(iii) In the Lagrangian formalism, explicit covariance is instead recovered. The first 
constraints one actually evaluates correspond to the secondary first-class constraints of the 
Hamiltonian formalism. At least at a classical level, the Lagrangian theory of constrained 
systems is by now a rich branch of modern mathematical physics [4,20,21], although the 
majority of general relativists are more familiar with the Hamiltonian framework. 

(iv) In the ADM formalism [5,6], the invariance group is not the whole diffeomorphism 
group, but a subgroup given by the Cartesian product of diffeomorphisms on the real line 
with the diffeomorphism group on spacelike three- surfaces. By contrast, in the Lagrangian 
approach, the invariance group of the theory is the full diffeomorphism group of four- 
dimensional Lorentzian space-time. This is the unique group responsible for the occurrence 
of constraints in Einstein's general relativity, on looking at symmetry properties of the 
action functional under such a group. 

(v) Jet-bundles theory provides a rigorous geometric framework for the hyperbolic 
problems of classical field theory and, in particular, general relativity. The corresponding 
evaluation of constraints, on using tetrad formalism, is elegant, very powerful, and well- 
suited for any attempt to study general relativity as a field theory with constraints. 

(vi) Although the elliptic boundary- value problems of Riemannian geometry and quan- 
tum gravity via Wick-rotated path integrals enable one to get a better understanding of 
different approaches to the quantization of gauge fields and gravitation [5] , the correspond- 
ing perturbative theory is non-renormalizable. Hence the background-field method should 
be complemented or replaced by a radically different view of space-time theories at the 
Planck length [22]. The Hamiltonian approach is the first step of the non-perturbative 
quantization program, but unfortunately it breaks covariance. 

We have thus tried to develop a classical multisymplectic analysis of Palatini and self- 
dual gravity theories which might be applied to a non-perturbative formalism for quantum 
gravity that preserves covariance. For this purpose, sect. 2 presents a derivation of Ein- 
stein's equations in tetrad form which relies on multisymplectic-geometry techniques (cf. 
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[23]). Section 3 studies multimomenta, sect. 4 derives constraint equations in multisym- 
plectic formalism, and their reduction to Ashtekar's form is obtained in sect. 5. The 
preservation of constraints is analyzed in sect. 6. Self-dual gravity is studied in sect. 7. 
Concluding remarks and open problems are presented in sect. 8. Relevant mathematical 
background is described in the appendix. 

2. - Multisymplectic form of Einstein's equations. 

Our analysis begins by studying the Palatini action Sp of (Lorentzian) general rela- 
tivity. This is a real- valued functional of the tetrad ej and the connection one- form ui a IJ , 
taking values in the Lorentz Lie-algebra, and given by 

Sp = \f d*xe e a j e b j nJ J . (2.1) 

With our notation, e = y/^g is the square root of the determinant of the space-time 
four-metric, and 

Q ab IJ ee d a u b IJ - d h u a IJ + uj K u b KJ - ^ K u a KJ (2.2) 

is the curvature of the four-dimensional connection one-form u a IJ . Moreover, a, b are 
tangent-space indices, whereas /, J are the usual internal indices [13]. 

In first-order theory, tetrad and connection are regarded as independent variables. 
Since we are aiming to use a Lagrangian version of first-order theory in terms of one-jet 
bundles (see appendix), it is useful to bear in mind that in general relativity one takes a 
fibre bundle whose base is space-time, and whose fibres are isomorphic to the Cartesian 
product of the space of Lorentzian four-metrics with the space of linear connections. In the 
language of tetrads and connection one-forms used in (2.1)-(2.2), one takes a fibre bundle 
Y which, in local coordinates, reads 
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To obtain the corresponding one-jet bundle J 1 (l r ) (see appendix), one is thus led to con- 
sider the mult i velocities V a bI corresponding to the tetrad, and the mult i velocities W ab IJ 
corresponding to the connection one- form. In local coordinates, our one-jet bundle J 1 {Y) 
is therefore represented by 



x a ,e a j,u a IJ ,V a bI ,W ( 



ab 



This leads to the Lagrangian 



L = g e ~i e J 



" " (W ab lj -W ba IJ + [u a ,u b ] IJ ) , 



(2.3) 



and hence to the Cartan four-form on J 1 (1 A ) as 

dL 



6r = L 



W u 



I J ^4 



<Px + 



dL 



du a IJ A d 6 x b 



e i e j - ej ej 



du a IJ A d 3 x b + uj K u b KJ d 4 x 



(2.4) 



Note that V a bI does not contribute, since derivatives of the tetrad do not occur in this 
first-order theory for vanishing torsion. The corresponding multisymplectic five-form fix, 
is obtained by exterior differentiation of ©x, as fix, = d<3>L- To write down field equations, 
one now considers a vector field U tangent to J 1 (1 A ). It has the form 



U = U° 



d 

dx* 



+ Uj 



d 



u 



d 



+ TJ d 8 + IT IJ () 



bl 



dW dl IJ • 



(2.5) 



One then takes the contraction of fix, with U. The pull-back of the resulting geometric 
object by means of the tangent lift j^(ip) of sections of the fibre bundle Y, leads to the 
Euler-Lagrange field equations. In our case one obtains 



iuVL L = -| (e a ! e b j - e b j e a j) e L h de\ A dw a IJ A d 2 x bm U r 



e 

+ 2 



(de a j)e b j + e",(de'j) - (de',)e'j - e',{de"j) 



b \„a 



A du a IJ A d z x b i U l 



+ 



( e - e \ - e \ e a j)e L h (de h L A d 3 x m U^ooJ K u t 



KJ 
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- \ [(de a j)e b j + e a Me h j) - (de b j)e a j - e b j{de a j)] uj K u b KJ A d 3 x, U l 
~ l(e a i e b j - e b j e a j) [{dwj K )u b KJ + uj K (dw b KJ )] Ad 3 Xl U l 



e 

+ 2 



U* e b j + U b j e a j - U b j e a ; - U a j e b j] 



du a IJ A d s x b + uj K u b KJ d 4 x 



-l^ a i^ b j-^ b i^j)e L h U h L 



e t e 



du a IJ A d 3 x b + uj K iv b KJ d 4 x 

a j)u a IJ e L h de h L Ad 3 x b 
(de a j)e b j + e a j(de b j) - (rfe 6 7 )e a 7 - e^cte",)] A d 3 x b 
+ \ {e a i e b j - e\ e a j) (uj K u b KJ + uJ K U b KJ )d'x . 



--u a IJ 

2 a 



(2.6) 



Such a lengthy equation is indeed necessary, since it involves many contributions which 
cannot be obvious to non-expert readers, and their interpretation needs a careful thinking. 
When one evaluates the pull-back of %u^l by means of the tangent lift 



(i (1) (^)) =(x\e%{x),uJ\x) d -^{x) dU « 



u 



dx b 



-(*) , 



the terms de h L and du a IJ occurring in (2.6) take the forms 



de\ = e h L j dx f 



du a IJ = u a IJ r dx l 



Thus, by using the identities 



dx f A d 3 xi = 6{ d 4 



x 



and defining 



dx f A dx h A d 2 x bc =(sl 5% - 5 f c 5^d 4 x , 



J) P h = TT m P h — TT h 

oe L — u e L m u L , 



6u a IJ = U™ uJ\ m -U a IJ , 



(2.7) 
(2.8) 

(2.9) 
(2.10) 

(2.11) 
(2.12) 
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L „a 



the field equations 



are found to take the form 



A>e a 7 = e% + u, b ^ e a L , 



(j«(^))*(^fi L )=0 



•(5e h L )G L h +{Sco a IJ )F a IJ = 0, 



where 



G\ ee 2e b j 



bh 



IL 



e t e 



d „L n U 



J ^ h bd 



F a u = D b 



e(e 6 / e a j - e\ e h j) 



Since e h L and u a IJ are independent, Eq. (2.15) implies 



(2.13) 



(2.14) 



(2.15) 



(2.16) 
(2.17) 



bh 



IL 



d L Q I J 
t J " h bd 



o, 



D h 



'(e a i 



i e j e i e j 



)] 







(2.18) 



(2.19) 



Eqs. (2.18) are the Einstein equations, while eqs. (2.19) express a property of a connection 
which is completely determined by the tetrad [13]. 

Note that in (2.11)-(2.12) the terms involving partial derivatives of the field variables 
are the horizontal part of the variation, and remaining terms represent the vertical part of 
the variation. Hence (2.11)-(2.12) have a very clear geometric meaning, and they lead to a 
considerable simplification of a lengthy calculation. Moreover, (2.13) defines the covariant 
derivatives of the tetrad. The connection D is a Lorentz connection which annihilates 
the Minkowskian metric 7]u on the internal space. The definition (2.13) is an additional 
condition we are imposing, since the action of D on space-time indices is not defined a 
priori [13]. 
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3. - Multimomenta. 



The analysis of constraint equations in sect. 4 makes it necessary to describe some 
properties of the multimomenta corresponding to the multivelocities defined in sect. 2. 
The multimomenta of a field theory are defined as the derivatives of the Lagrangian with 
respect to the multivelocities. In general relativity, the multimomenta resulting from the 
Lagrangian (2.3) are defined as the densities 

81 

H» S i=C, (3.2) 

UV bl 

where the p ab jj are bivectors defined as 

P ab u = e% e h j - e\ e a j . (3.3) 

Of course, the multiplicative factor in (3.1)-(3.2) is unessential, and is introduced for con- 
venience. Note that the multimomenta ir a bI vanish, and this reflects that torsion vanishes 
in general relativity. 

From the definition (3.3), we note that the bivectors p ab u satisfy the following com- 
mutation relations: 

"p.* p cd] ^ = p ac u gbd + ^ g ac _ pad j g bc _ p bc ^ j g ad (3 ^ 

In particular, on a spacelike hypersurface, Eq. (3.4) may be used to derive the identity 

^uf-A^-AjP' . (3-5) 



p l0 ,p>° 



J ij 



where the index refers to the time coordinate and the indices i,j refer to the spatial 
coordinates. This property will be useful in sect. 5. The multimomenta correspond to 36 
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variables (6 for each bivector density), therefore equation (3.5) expresses the eighteen vari- 
ables associated to the V 3 u as quadratic functions of the p 0z jj bivector densities. In the 
following sections the constraint equations are all expressed in terms of the multimomenta. 



4. - Constraints. 



When one studies classical mechanics and classical field theory one learns that, if the 
Lagrangian is invariant under the action of a group, then by virtue of Noether's theorem 
there exist functions which are constant along solutions of the equations of motion. In the 
case of the invariance under a gauge group or the diffeomorphism group of general relativity, 
such first integrals always vanish along solutions of the equations of motion. Hence the 
first-class constraints of a field theory result from Noether's theorem through the action of 
the gauge group or the group of space-time diffeomorphisms [23,24]. The multimomentum 
map is the mathematical tool which enables one to describe these properties of classical 
fields. With the notation of Eqs. (A.1)-(A.2), the multimomentum map on a section of 
our jet-bundle J X {Y) is defined by the expression [23,24] 



(i (1) (¥0)*[^r(j c ) 



dL 



d 3 x f 



(4.1) 



where is the variation along the fibre and £ a is the variation along the base space. In 
terms of the variables e a 7 and u a IJ the multimomentum map J becomes 



m = 



dL 

'de a T 



IJ C c 



d"a IJ ,c 



d 3 x c . (4.2) 



In the Hamiltonian framework, setting to zero the integral of the multimomentum map 
on a spacelike hypersurface E leads to the first-class constraints of the theory [23,24]. 
Moreover, in the case of null hypersurfaces, second-class constraints also occur [25], and 
hence they cannot be derived from the multimomentum map. On passing from spacelike 
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to null hypersurfaces, discontinuities occur in the normal to the hypersurface and in the 
induced three-metric [25]. Hence there might be a covariant constraint analysis for any 
spacelike hypersurface, and an independent, covariant constraint analysis for any null 
hypersurface. [We are grateful to Luca Lusanna for bringing this open problem to our 
attention] 

To express our J(£), from 

dL 



r)p a 
ue I,c 



0, 



and 



L U = ~t\a + (D a X) IJ , 

it follows (on using (3.3)) 

Wt] = [ AO = \\ \ep ac u{e, a " b U - (Da A) IJ W,6 ?) 



(4.3) 



(4.4) 



(4.5) 



where one has defined the covariant derivative with respect to the tetrad indices (cf. (2.13)) 



(D a X) IJ = (d a \) 1J + [\,u a 



u 



(4.6) 



Thus, taking the integral of the multimomentum map on the spacelike hypersurface E, 
integrating by parts for the term in A/ 7 , and defining 



ac — ~ac \IJ 

a = p u a 



(4.7) 



one obtains 



= \J X IJ (D a p ac )u d 3 x c + l -J^ 



d n a ac d s x r 



(4.8) 



Again from integration by parts in the second integral, and defining 

ac — ~ac IJ th 

P = P I J U b S ) 

10 
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one finds 



Ht] = \\ -i b (p aC Ij) (V J ,a " "a U , b ) ~ e P aC IJ,a ^ + b^IJ *lJ J £ 

J X] l 

+ \ X + xU ( D aP ac )iJ d 3 x c d a a ac d 3 x c . 



d 6 Xr 



(4.10) 



Using equation (2.19), and denoting by iff s [£] the part of /e[£] not involving total diver- 
gences, one can show that 

1 



£ d d 3 x c 



^j^\ IJ {D a p ac )jj d 3 x c -^J^Tr 



P ac tt ad - \p ab Sl ah 5 c d 



£ d d 3 x c , (4.11) 



where X IJ and £ d are independent and arbitrary quantities. Then from imposing /e[£] = 0, 
one finds two sets of constraints: 



/ d a a ac d 3 x c - [ X IJ (D a p ac )ud 3 x c = 0, 



(4.12) 



and 



d a p ac d 3 x c - / Tr 



1 

r 



P ac n ad - ^p ab n ab 5 c d 



i d d 3 x c = 



(4.13) 



Note that the total divergences appearing in (4. 12)- (4. 13) lead to boundary terms evaluated 
on the two-surface 9E. Hence sufficient conditions for the vanishing of such boundary terms 
are as follows: (i) E has no boundary; (ii) p ab jj vanishes at <9£; (iii) X IJ vanishes at <9E (in 
(4.12)); (iv) £ b or u b IJ vanishes at <9E (in (4.13)). Although these conditions are (rather) 
restrictive, from now on we will always assume that (i) or (ii) or both (iii) and (iv) hold. 
Hence (4. 12)- (4. 13) reduce to 



/ X IJ (D a p ac ) 
Jt, 



ij d x c = , 



(4.14) 



and 



L 



Tr 



p aC ttad - \p ab Sl ab S d 



^ d 3 x c = / G L d e\ i d d 3 x c = . 



(4.15) 
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5. - Reduction of constraints to Ashtekar's form. 



In this section we show that the previous equations reproduce Ashtekar's results for a 
Palatini Lagrangian [13]. For this purpose, we choose the adapted local coordinates defined 
by the condition x° = constant, while the remaining three coordinates are denoted by the 
spatial indices i,j, k. In view of the covariance of the Lagrangian formalism, it is enough 
to prove the equivalence with Ashtekar's constraint analysis on spacelike hypersurfaces in 
this particular coordinate system. The constraints (4. 14)- (4. 15) are then found to take the 
form (from the arbitrariness of X IJ and £ d ) 



(D a p a0 )u = 



(5.1) 



which corresponds to the Gauss constraint [9,13], and 



Tr 



P a0 n ad - ±p ab n ab s° d 



= o 



(5.2) 



Note that (5.1) reflects the invariance of general relativity under local Lorentz transforma- 
tions, and leads to six independent constraint equations. They may be further split into 
internal rotations and boosts [9,13]. The last four constraints (5.2) can be split into two 
parts, for d ^ one has 



Tr 



p a0 n 



aj 



= Tr 



p i0 n 



= o, 



(5.3) 



which can be identified with the vector constraint; and for d = 



Tr 



/7' l, a„, - \p ab n ab 



Tr 



= Tr 

2 



1-3 







(5.4) 



which is the Hamiltonian constraint. To see that (5.4) implies the Hamiltonian constraint 
in its usual form, we point out that from (3.5) 



Tr 



F'SUj =2e-\g 00 )- 1 Tt p i0 fP Q tJ +2(^ 00 )- 1 ^°Tr 



(5.5) 
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Then from (5.3)-(5.5) it follows that 



Tr 



?V°% I = o , 



(5.6) 



which is the familiar constraint equation quadratic in the momenta [13]. 



6. - Preservation of constraints. 



We have now to prove explicitly that the constraints (4.14)-(4.15) are preserved in our 
multisymplectic approach. For this purpose, we begin by requiring that the following Lie 
derivative: 



1 = 



= j l v Tr(p ac n ad -U ab n ab s^jc 



d 3 x c , (6.1) 



should vanish, where L v denotes the Lie derivative along a smooth vector field r\. The 
integral (6.1) results from taking differences of the constraints evaluated at two generic 
spacelike hypersurfaces E and related by a one-parameter flow generated by rj. Thus, 
since the Lie derivative along 1] of a weight-w vector density V reads [26] 

(L„y) ° = d b (V a v b - FY ) + V a d b V b + (w- l)V a d bV b , (6.2) 

one finds that (6.1) takes the form 

1 = 2 f G c d £V d 2 x bc + [ rfG b c V 6 C C d 3 x a + f rf f V b G b c d 3 x a , (6.3) 

where G a h = e G L b e a L (see (2.16)), and V is the (torsion-free) connection which annihi- 
lates the tetrad. Hence V a V a = d a V a [26]. The first term in (6.3) vanishes by virtue of the 
boundary conditions imposed in sect. 4, and the second integral vanishes along solutions 
of the field equations. Last, the third integral vanishes by virtue of the contracted Bianchi 
identities with respect to the space-time connection V. 



13 



SPACE-TIME COVARIANT FORM OF ASHTEKAR'S CONSTRAINTS 

Similarly, one finds that also the Gauss-law constraint (4.14) is preserved, since 

J L v X IJ (D a p aA ) ^ d 3 x c = 2 J V b X IJ (D a p ac ) ^ d 2 x bc 

+ J rf (<9 b A 7 J ) ( D a p a6 ) ^ d 3 x c + J rf X IJ d b ( L> ai 5 a6 ) ^ d 3 x c . (6.4) 

Again, the boundary term vanishes by virtue of the boundary conditions of sect. 4, and 
the second term vanishes on imposing the field equations. Moreover, the third term on the 
right-hand side of (6.4) vanishes by virtue of the identity 

d b (D a p ab ) =l\n ba ,p ab ] =o. (6.5) 

V / u 2 L J i j 

Another way to derive these results is to consider the four-dimensional volume inte- 
gral of V 1 c {g c d taken over a region V whose boundary is the disjoint union of two 

hypersurfaces E and £'. On applying the Leibniz rule and imposing the field equations, 
one finds the equation 

/ G c d e d d 3 x c - [ G c d e d d 3 x c = I t d V c G c d d 4 x . (6.6) 

Again, the preservation of constraints is achieved by virtue of the contracted Bianchi 
identities. A similar argument can be applied to prove the preservation of (4.14) (cf. 
(6.4)). 

In a Palatini formalism, one has also to consider the extra constraints corresponding 
to second-class constraints [13]. In the Hamiltonian formulation, second-class constraints 
arise since one starts from tetrads and connection one-forms, while the constraint analysis 
makes it convenient to replace the tetrads by the momenta. To prove equivalence between 
these two formulations it is then necessary to restrict the momenta so as to recover general 
relativity. One then finds six primary second-class constraints and six secondary second- 
class constraints [13]. This leads to the two degrees of freedom of real general relativity 
(page 59 of ref. [13]). A similar problem occurs in our analysis, if the action (2.1) is re- 
expressed in terms of the multimomenta. One has then to consider additional conditions 
obeyed by the multimomenta. 
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In the light of (3.1) and (3.3), the second-class constraint equations found in ref. [13] 
take the form 

n c n d e 1 JKL p ac LJ p b d KL = , (6.7) 



on on J JKL ~cd M 

n d nf e p j 



P af M j D ^hP bh ) KL +P bf M j D c (n h p ah ) 



KL 



, (6.8) 



where n a is the unit timelike normal to the spacelike hypersurface E. Note that our 
constraint equations (4.14)-(4.15) and (6.7)-(6.8) are all expressed in terms of the multi- 
momenta. 

At a Lagrangian level, however, there are no primary constraints [4,20], since one 
deals with the pull-back on a manifold corresponding to the primary-constraint submani- 
fold of the Hamiltonian formalism. The problem remains to derive the constraints (if any) 
corresponding to the secondary second-class constraints of the Hamiltonian formalism (of 
course, constraints cannot be divided into first- and second-class in a Lagrangian frame- 
work, since no Poisson brackets exist) . They cannot be derived from the multimomentum 
map, which only applies to the analysis of constraints which are a counterpart of first-class 
constraints. In ref. [20], we were able to use the Gotay-Nester Lagrangian analysis in a 
finite-dimensional model to derive the Lagrangian counterpart of secondary second-class 
constraints. However, we do not yet know whether that technique can be extended to our 
formulation of general relativity. 

For this purpose, we are currently investigating a more general set of equations obeyed 
by the multimomenta, i.e. (cf. [26,27]) 

tabcd P ab u p Cd KL = T[UKL] , (6.9) 

where T is a tensor density proportional to €jjkl- Eq. (6.9) admits, as a particular case, 
the condition of simplicity of the multimomenta (e.g. set / = K, J = L in (6.9)) in the 
abstract indices a, 6, as well as the analogous of (6.7) where the roles of space-time and 
internal indices are interchanged, i.e. n J n L T[ IJKL ] = 0, where n J = e a J n a . 
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7. - Self-dual gravity. 

The Hamiltonian formulation of self-dual gravity has received careful consideration 
in the current literature [13]. The corresponding formalism, however, is not manifestly 
covariant. To overcome this problem, Samuel [28], and, independently, Jacobson and 
Smolin [29,30], proposed a Lagrangian approach based on a self-dual action. The key idea 
is to take complex (co)tetrads on a real Lorentzian four-manifold [31], and then express 
the action functional in terms of the tetrad and of the self-dual part of the connection: 




(7.1) 



Remarkably, the complex self-dual action 

S SD = \ I d 4 x e e aI e bJ n abIJ (+u) (7.2) 

is related to the real Palatini action (2.1) by [32] 

S SD [e, + uj(e)} = ls P [e,u(e)] - % - f d 4 x e e aI e T JKL Vt aJKL {u{e)) , (7.3) 
1 ° J M 

where the second term in (7.3) vanishes for vanishing torsion, since then Q a [jKL] = 0. 
Hence the resulting field equations for real general relativity are equivalent, and the corre- 
sponding constraints are first-class only (before imposing reality conditions) and polyno- 
mial (page 64 of ref. [13]). 

The self-dual equations are hence obtained by replacing the full connection with the 
self-dual connection in sees. 2, 4, 5 and 6. Thus, on defining 

+ a ac = +pa Cij X IJ ^ ( 7 4 ) 

+ p ac = + p ac u + u; b IJ e , (7.5) 
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the constraint equations become (cf. (4. 12)- (4. 13)) 

/ d a +a ac d 3 x c - [ \ IJ (D a + p ac ) u d 3 x c = , (7.6) 
Jt Jt 

[ d a +p ac d 3 x c - [ Tr 
Jt, Jt, 

In other words, in the most recent presentations [13,31], the equivalence between general 
relativity with a Palatini action and its self-dual version is proved. In this paper we 
follow the same argument and we find entirely analogous results, while the multisymplectic 
framework enables one to preserve covariance. Moreover, no extra constraints are necessary 
to recover the original content of the theory, when the action (7.2) is re-expressed in terms of 
the multimomenta (cf. [13]). As a last step, one has to impose suitable reality conditions to 
recover real general relativity [13,31]. Their formulation in the multisymplectic framework 
is not studied in our paper, and is a subject for further research. 

We think one should emphasize again that, in the canonical-gravity approach to self- 
duality, one takes complex (co)tetrads on real space-time four-manifolds. By contrast, 
in other branches of modern relativity [33], one is interested in four-complex-dimensional 
complex-Riemannian manifolds. In such a case, no complex conjugation can be defined, 
since this map is not invariant under holomorphic coordinate transformations, and no 
four-real-dimensional sub- manifold can in general be singled out [33]. Hence the problem 
of reality conditions to recover real general relativity cannot even be addressed in the 
complex-Riemannian framework. The corresponding theory of self-duality involves the 
Weyl spinors, and is not equivalent to the model outlined in this section [33]. 

8. - Results and open problems. 

This paper has studied Palatini and self-dual gravity theories by using tetrad formal- 
ism and multisymplectic techniques in Lorentzian four-manifolds (cf. [34]). Our results 
are as follows. 



+p ac n ad (+u) - ±+p ab n ab (+u) s c d 



£ d d 3 x c = 



(7.7) 
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(i) The first-class constraint equations of Palatini theory are given by (4.12)-(4.13). Inter- 
estingly, boundary terms occur, and they vanish under the sufficient conditions listed at 
the end of sect. 4. 

(ii) The Hamiltonian constraint of general relativity is linear in the multimomenta, as 
shown in (5.3)-(5.5). Indeed, on studying Palatini formalism, Eq. (3) of chapter 4 of ref. 
[13] implicitly expresses this property. Our analysis, however, makes it more evident. The 
multisymplectic framework of sect. 2, and the multimomenta formalism of the following 
sections, seem to add evidence in favour of the Lagrangian point of view being able to 
supplement the Hamiltonian formalism to get a better understanding of the gravitational 
field. Moreover, the linearity in the multimomenta of all constraint equations may have 
far reaching consequences for the quantization program. 

It should be emphasized that our analysis has not improved the understanding of the 
initial-value problem in general relativity. However, our equations are covariant in that 
they do not depend on a particular spacelike or null hypersurface, and they naturally lead to 
study the space of multimomenta. Hence they are a step towards a covariant formulation of 
relativistic theories of gravitation regarded as constrained systems. Within this framework, 
it may be interesting to analyze 2+1 gravity and theories with no n- vanishing torsion. 

A naturally occurring question is whether our classical analysis can be used to formu- 
late an approach to non-perturbative quantum gravity which improves the results obtained 
within the more familiar Hamiltonian framework [13]. For this purpose, it appears neces- 
sary to get a better understanding of covariant Poisson brackets [35], and possibly of the 
geometric quantization program [36,37]. Hence we do not yet know the key features of 
the resulting quantum theory. However, the elegance of the mathematical formalism, and 
its wide range of applications at the classical level, make us feel that new perspectives in 
non-perturbative quantum gravity are in sight. 
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Appendix. 

Since our paper is primarily addressed to physicists interested in general relativity, we 
limit ourselves to a very brief outline of some geometric ideas used in our investigation. 
An extended treatment may be found in ref. [23]. 

In our paper, the notation J 1 (Y) means what follows [24]. Let A be a manifold and 
let Y be a fibre bundle having A as its base space, with projection map ttxy ■ A fibre is 
then given by n^} Y {x), where x G X. Moreover, let 7 : T X X — > T y Y be a linear map 
between the tangent space to A at x and the tangent space to Y at y G 7r^y(x). Now, 
given a point y belonging to the fibre Y x through x G A, we consider all 7 maps relative 
to y G Y x . This leads to a fibre bundle J X {Y) having the fibre bundle Y as its base 
space and fibres given by the 7 maps. Such a J 1 (Y) is called the one-jet bundle on Y. 
If Lp( A \x^) is a section of Y, the tangent lift of ip^ to a section of J l {Y) is denoted by 
j^ 1 ' (</?). It is given by the map 

^>:^(^V),*^). (Al) 
The Legendre map, whose pull-back appears in Eq. (4.1), is a function 

TL : J\Y) ->|>(y)]* , (A2) 



where 



J 1 (y ) is a fibre bundle having Y as its base space, and whose fibre through y <EY 



is given by the affine maps on the elements of the fibre of J 1 (^), with coefficients in the 



bundle of n + 1 forms on 1 at 1. Hence 



is called the dual of the one-jet bundle 



J 1 00 

J 1 (y). If, in local coordinates, J X {Y) is described by ^ , the expression of its 

dual in local coordinates is given by y^ A \p,p^ A ^ , where the definitions of momenta 
and of Legendre map yield 
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p^L-p^v^. (AA) 

* * * 
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